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1. The mathematical theory of the deformation of a thin elastic plate or shell 
involves difficulties of a formidable nature. This is partly owing to the fact that 
an approximate solution, which is correct as far as terms involving the cube of 
the thickness of the plate or shell, cannot be obtained without having recourse 
to an extremely lengthy and troublesome process which requires the Calculus of 
Variations, and partly because the older writers upon this subject based their 
investigations upon hypotheses which were in most cases inadequate and erro- 
neous. "When a mathematician of standing and reputation gravely propounds a 
hypothesis which turns out to be incorrect, or condemns as unsound some 
method which is perfectly legitimate, and in addition obtains by means of his 
erroneous hypothesis many results which are substantially correct, the mischief 
done to science can hardly be exaggerated; for subsequent investigators not only 
are led astray from the path which must be followed in order to obtain a satis- 
factory theory, but are hampered by the difficulty of convincing mankind of the 
errors of their predecessors. 

The development of the theory of thin elastic plates and shells has been 
retarded by two errors for which Olebsch and Saint-Venant are mainly respon- 
sible. The first error is that the three stresses B, S, T are accurately zero 
throughout the substance of the plate or shell ; the second error is that it is not 
permissible to expand the various quantities involved, in a series of ascending 
powers of the distance of a point from the middle surface. In England the first 
error may, I think, now be regarded as exploded ; but the second one has not 
yet been completely driven from the field, for a perusal of Mr. Love's recent 
Treatise on Elasticity shows that he entertains some bias against the method 
of expansion, though on what grounds I am at a loss to conceive. 

Owing to the unsatisfactory manner in which Mr. Love has dealt with the 
theories of thin plates, shells and wires in his book, I propose in the present 
paper to give as concise an account of the first two theories as the nature 
34 
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of the case will admit. I shall commence by expounding the fundamental 
principles of these theories ; I shall then proceed to develop them in a form 
suitable for mathematical calculation ; I shall endeavor to avoid all unnecessary 
mathematical complications; and I shall omit or pass over very lightly those parts 
which are not of much physical interest, or which require lengthy and difficult 
analysis for their elucidation. 

2. We shall first consider the resultant stresses which act upon any element 
of a thin shell. 

C 




Let OADB be a small curvilinear rectangle described on the middle surface 
of a thin shell, the sides of which are lines of curvature ; and let us consider a 
small element of the shell bounded by the external surfaces, and the four planes 
passing through the sides of this rectangle which are perpendicular to the middle 
surface. 

The resultant stresses per unit of length which act upon the element, and 
which are due to the action of contiguous portions of the shell, are completely 
specified by the following ten quantities, viz. across the section AD : 

7i = a tension across AD parallel to OA, 

M % = a tangential shearing stress along AD, 

iV a = a normal shearing stress parallel to OG , 

G % = a flexural couple from G to A whose axis is parallel to AD, 

H x = a torsional couple from B to G whose axis is parallel to OA . 

Similarly the resultant stresses per unit of length which act across the section 
BD are : 

T % = a tension across BD parallel to OB, 

M 1 = a tangential shearing stress along BD, 

Ni = a normal shearing stress parallel to OG, 

6?i = a flexural couple from B to whose axis is parallel to BD , 

H % = a torsional couple from G to A whose axis is parallel to OB. 
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By resolving all the forces and stresses which act upon the element parallel 
to OA, OB and 00, and by taking moments about these lines, we obtain the six 
equations of motion of the element ; but these equations will not enable us to 
solve any statical or dynamical problems, since they furnish only six equations 
connecting ten unknown quantities. In order to complete the solution we 
require the values of the stresses and the equations of motion in terms of the 
displacements of a point on the middle surface and their space variations in 
terms of the coordinates of that point. 

The internal stresses at any point of the substance of the shell are com- 
pletely specified by the six quantities P, Q, B, S, T, U; and at this stage of 
the investigation it will be well to state that the notation of Thomson and Tait's 
Natural Philosophy will be employed for internal stresses and elastic constants, 
but that the three extensional strains parallel to OA, OB, 00 will be denoted 
by Gu 2 > ff 3> an d the three shearing strains about these lines by 5 1( 5 2 , o s . 

3. In most cases of practical importance (although there are a few excep- 
tions), the surfaces of the shell are free from stress, and when this is the case, 
the three stresses B, S and T must vanish at both surfaces. Now Clebsch and 
several of the older writers fell into the mistake of supposing that these stresses 
are zero throughout the substance of the shell, and by the aid of this erroneous 
hypothesis theories were constructed which nevertheless furnished numerous 
results which were substantially correct. The reason of this was as follows : 
The theory of thin plates and shells is an approximate one, in which the thick- 
ness is supposed to be a small quantity in comparison with the least radius of 
curvature of the middle surface, and it is sufficient to carry the approximation 
no farther than powers involving the cube of the thickness. If 2h is the thick- 
ness and hi the distance of any point of the shell from the middle surface, it will 
appear as we proceed that if the lowest terms in B , S, T were quadratic func- 
tions of h and h', the retention of these stresses would give rise to terms of a 
higher order than h s , and consequently an approximate solution could be 
obtained to that order by treating B, S, T as zero. The first question therefore 
is, have we any evidence which would justify us in concluding that this supposi- 
tion is true ? 

4. The periods of the vibrations of a plane plate of infinite extent have 
been worked out by Lord Rayleigh* by means of the general equations of an 

* Proc. Lond. Math. Soc, Vol. XX,.p. 225. 
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elastic solid. The vibrations consist of two distinct types, viz. flexural vibra- 
tions unaccompanied by extension of the middle surface in which the displace- 
ment is perpendicular to that surface, and extensional vibrations in which the 
displacements are parallel to the middle surface. If Lord Rayleigh's expres- 
sions for these two kinds of vibrations be expanded in powers of the thickness, 
it will be found that if Injpi, ZnlPi be the periods, 

^ = 3 P (T+ro) + — higher P owers of A2 > (l ) 

where n is the rigidity, m = h -\-\n, where h is the resistance to compression, 
p the density and 27i//the wave-length of the disturbance. The value of^? g is 

±mnf f {m-nffh? ) 

Pi ~ p(m + n) \ 1— 3 (m + nf + ' ' * ' j ^> 

Now if we attempt to obtain the same results by assuming that the three 
stresses R, S, T &xq quadratic functions of h and h', it will be found that these 
quantities give rise to terms in the expression for the energy which are of a 
higher order than h s and may therefore be neglected, and we shall finally obtain 
the same expressions for p x and p % as are given above.* 

It can also be proved by means of Lord Rayleigh's paper, that in these two 
particular cases R , S and T are quadratic functions of h and h 1 . 

In order to test the truth of the hypothesis in the case of a curved shell, I 
workedf out the period of the radial vibrations of a complete cylindrical shell to 
a second approximation, so as to obtain the term in h 3 , by means of the general 
equations of elasticity and also by means of the theory of thin plates, and I 
found that the results obtained by these two methods agreed. The value of p 
in this case is 

g> - 4wm f A* r 4(iw -n) -i ] () 

where a is the radius of the middle surface. 

We notice in passing that in the case of a cylindrical shell the pitch of the 
purely extensional vibrations rises as the thickness increases, whereas the con- 
trary is the case when the surface is an infinite plane plate. Also the frequency 
of the flexural vibrations of a plane plate is proportional' to the thickness, 
whereas the frequency of the extensional vibrations contains a term independent 

*See Proc. Lond. Math. Soc., Vol. XXI, p. 51. t Ibid. p. 53. 
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of the thickness ; the notes yielded by extensional vibrations are therefore of 
much higher pitch than those yielded by non-extensional ones. 

5. In consequence of the direct evidence furnished by the above-mentioned 
results, I felt justified, in 1889, in laying down the following fundamental 
hypothesis as the basis of the true theory of thin plates and shells, viz. : 

When the surfaces of a thin plate or shell are not subjected to any surface 
forces, such as external pressures or tangential stresses, the three stresses R, S, T, so 
far as they depend upon h and h', are of the form 

ih + u 3 "H • •••«»+••••, 

where u n is a homogeneous n ti0 function of h and h'. 

On applying this hypothesis to obtain a theory of plane plates and cylindri- 
cal and spherical shells,* I found that expressions for the kinetic and potential 
energies could be obtained which are correct as far as h 3 , and that in deducing 
these expressions a knowledge of B, S and T was not required, but that they 
might be treated as zero, inasmuch as they would, if retained, give rise to terms 
involving powers of h higher than h s . I was also able to test the correctness of 
the fundamental hypothesis and also of the work in the following manner: — 

As soon as the correct expressions for the potential and kinetic energies 
have been obtained in terms of the displacements, the equations of motion and the 
values of the ten sectional stresses can be deduced in terms of these quantities 
by means of the Principle of Virtual Work. 

Now we have already pointed out that the six equations of motion can be 
written down in terms of the ten sectional stresses, and if we substitute in these 
equations the values of these stresses in terms of the displacements which we 
have obtained by means of the Principle of Virtual Work, three of them ought 
to be identically satisfied, whilst the remaining three ought to be identical with 
the three equations of motion which we have obtained by means of that prin- 
ciple. This is the first test, and it is satisfied. The second test is that all the sec- 
tional stresses except T x and T % can be calculated by a direct method, and the 
values so obtained ought to agree with those furnished by the Principle of Virtual 
Work. This test is also satisfied. 

6. Let F be any point of the shell, P its projection on the middle surface ; 
then all the quantities with which we are concerned are functions of the position 

*Phil. Trans., 1890, p. 483. 
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of P, and are therefore functions of (r , z, $) or (r, 6, ty) according as the shell 
is cylindrical or spherical. If, therefore, C^' be the value of any such quantity 
at P, and ^ its value at P, it follows that 

ty>=F(r) = F(a + h') 

-«+<£)+**■($)+ ■■•■ 

by Taylor's theorem, where the brackets are employed, as will be done through- 
out this paper, to denote the values of the differential coefficients at the middle 
surface. 

I have already alluded to the fact that objections have been raised by 
Saint- Venant, and endorsed by Mr. Love,* against the method of expansion. 
These objections appear to me to be so utterly without foundation that it 
would be scarcely worth while to labor the point were it not that Mr. Love, 
throughout his recent treatise, has shown a persistent disinclination to employ 
this method, although he does not venture to repeat his former arguments. In 
consequence of this he has been compelled to resort to lengthy and unnecessarily 
complicated methods which I fear will greatly retard the further development of 
the subject. It has been known for many years past that when a function and 
its differential coefficients are finite and continuous between given limits of the 
variable, the function can always be expanded between these limits by Taylor's 
theorem, and that we can stop at any term of the series that we please and 
express the remainder in a finite form. Now although it is quite true that 
functions may be imagined which, by reason of their becoming infinite or dis- 
continuous between the limits, would render the application of Taylor's theorem 
inadmissible, yet the occurrence of such functions is precluded by physical con- 
siderations, because they would involve a rupture of the material. We all know 
that bells sometimes crack and wires break ; but such catastrophes are not con- 
templated, and would require a special theory to account for them. The sooner 
these objections are consigned to oblivion the better it will be for the theory of 
Elasticity. 

We are now in a position to consider the mathematical part of the theory. 

Theory of Plane Plates. 

7. We shall commence by obtaining the equations of motion in terms of 
the sectional stresses. 

•Phil. Trans., 1888, p. 491. 
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Let u', xf, w' be the displacements parallel to OA, OB, 00 of any point 
P' of the plate ; X, T, Z the components of the bodily forces per unit of mass. 
Then taking moments about 00, it follows at once that M x = i/ 2 ; also resolving 
parallel to OA, OB, 00 and taking moments about OA, OB, we obtain the 
following five equations, viz. : 



dT x dM 

dx ' dy 
dT . dM 



dy 

dK 



' dx 



+ 



dN 1 



dx dy 

dG x , dH t 



dy 



+ 



dx 



+ Ni =■ 



dG % dR Z jrjr _ 



dx 



dy 



/h 
vldh' — 2p7LZ, 

ofi'dh'—2ohY, 

pfw'dh'—2phZ, 

- pf&h'dh', 

pfit'h'dh'. 



(1) 



8. We shall now proceed to calculate the four couples in terms of the dis- 
placements of a point on the middle surface. 

Let P be the projection of P on the middle surface, and let the values of 
the different quantities at F and P be distinguished by accented and unaccented 
letters. The equations determining the six strains in terms of the displace- 
ments are 



%: 



du' 
dx 
dv' 



dz dy 



dtd 



dv 1 
'dy 
did 



dvtf 



<t; = 



° 3== ^ + 



du' 

dz 



o 3 






dz 

du' , dv' 



dy 



+ 



dx 



(2) 



The value of B', when expanded in powers of h', may be written 
B' — A + A 1 h' + $A 2 h'*+ 



(3) 



where, according to the fundamental hypothesis, A and A x do not contain any 
lower powers of h than h z and h respectively. But 

B' = (m + n) a' s + (m — n)(a[ + a' t ) 
= (m + n) <t 3 + (m — rc)(oi + <r,) 

+ \( m +n)(^f)+(rn-n)(Z, + ti)\h> +.-..., (4) 



(5) 
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where % and p denote the values of dajdz and da % jdz when s= 0. Accord 
ingly, by (3) and (4), 

(m +n)o 3 + (m — n)^ + ff 2 )= A , j 
(m + n)(**) + (m - n)(a + p) = 4, J 
The value of 6?! is 

ffi = —jy«M =-w <^§) , 
(f) =(-+-)».+(—) | *+&)}. 

whence putting 



*=^, (6) 

we obtain 

G £ 1 = -|»A 3 ] l « + J E'(^ + ^) + ^A/2«}. (7) 

Similarly, G 8 = fnft 8 {3, + ^(3, + ft) + EA 1 /2n\. (8) 

The value of A x is unknown, but according to the fundamental hypothesis 
A x does not contain any lower power of h than the first, whence the terms in A 
must be neglected, since the approximation is not carried farther than h 3 . 

Now by (2) 

„ da 1 d?u da 2 cPw 

dz dzdx dx dx 2 ' 

also _ dcSi d?v doi dPw 

P" dz dzdy dy dy % ' 

Since S' = na[ and T' = nu! 2 , it follows from the fundamental hypothesis that 
au a it which are the values of ai, g^ at the middle surface, do not contain any 
lower power of h than W ; accordingly the terms involving these quantities may 
be neglected when multiplied by h 3 , whence the values of O t , G % finally become 



where v 2 has its usual meaning. 
Again, 

«=-n=-i*(f) 



(») 
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and da s d?u d?v da% du^ d?w 

dz dydz dxdz dy dx dxdy ' 

so that to the above order of approximation 

B.=-n,=*»£%. do) 

Equations (9) and (10) completely determine the four couples. 

9. If we expand vl, v 1 , «/ in powers of h' and eliminate their differential 
coefficients with respect to z, the right-hand sides of the last three of (1) can be 
calculated by a similar process, and we shall obtain 



dy + dx + Ml - T P* dy ' 
dGn , dH % , T 3 -,-dw 

•35- + -^-^ =- i ? h d^- 



(10, A) 



Substituting the values of the couples from (9) and (10) in the last two, 
and the resulting values N t and N 2 in the first, we obtain the following equation 
of motion, viz. : 

inh i (l + E)V i w+ 2p(io — Z) + %ph 2 (E — 2)V% = 0. (11) 

If Z= the last term must be omitted, since it gives rise to a term A 4 in 
the period, and (11) becomes 

4^_ V ^+p^=0. (12) 

3 (m -f- n) r 

The motion represented by this equation consists of a displacement which 
is perpendicular to the middle surface ; moreover, the motion does not involve 
any extension or contraction of this surface, for the extensional strains are given 
by the first, second and sixth of (2), and when u and v are zero these strains are 
also zero. Vibrations of this kind are sometimes called lateral vibrations on the 
ground that the displacement is perpendicular to the middle surface, and some- 
times flexural vibrations,, because they involve flexion or bending of this surface 
unaccompanied by extension. The periods of these vibrations are much longer 
than those which depend upon extension alone, and they constitute the most 
35 
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important class of vibrations which a plane plate is capable of executing. They 
are fully discused in Vol. I, Chap. X, of Lord Rayleigh's Theory of Sound. 

10. We must now consider the boundary conditions which hold good at a 
free edge. 

It is obvious that the flexural couple about a free edge must vanish ; hence 
if the free edge is parallel to the axis of y, # 2 = 0, and therefore by the second 
of (9), 

It might also be thought that the torsional couple B and the normal shear- 
ing stress N must also vanish, in which case we should have three boundary 
conditions. This was Poisson's view, but Kirchhoff showed that only one addi- 
tional equation was necessary, and it is now well known that Poisson's boundary 
conditions are erroneous. 

The reduction of Poisson's boundary conditions depends upon the curious 
theorem that it is possible to apply a system of stresses to the edge of the plate 
without doing any work. By Stokes' theorem 

the integration extending round the rectangle OADB. 
If we apply to the side AD the stresses 



2r- dJsr 


H[ = B?, 


W-dS' 


m=- 



to the side DB the stresses 

JTTI 

3', 

and to the sides BO, OA corresponding and opposite stresses, the preceding 
integral becomes 

which shows that the work done by this system of stresses is zero. Such a 
system of stresses may therefore be applied or removed without interfering with 
the equilibrium or motion of the plate. 
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If, therefore, we suppose that the rectangle OADB, instead of being under 
the action of the remainder of the plate, is isolated, and that constraining stresses 
are applied to its edges, then it follows that if instead of the torsional couples 
.Hi, — H x (since 5^= — ZZi) due to the action of contiguous portions of the 
plate, we apply torsional couples fjjj, J| 2 , where 

1L= fli + -ff', 



--H 1 -ff', 



(13) 



the energy and state of strain will be unaltered, provided we apply in addition 
the stresses 

= * + ** 



dy ' 

AT _L dH ' 

i = ^ + WJ 



(14) 



whence eliminating H' between the first and the, second of (13) and (14) respec- 
tively, we obtain 

gUi _ xr dHt 

dy 

dl 



N t 



*+-% = * + 



dy ' 
dH, 



dx ' , 



(15) 



In these equations the Roman letters denote the values of the stresses due 
to contiguous portions of the plate, whilst the Old English letters denote the 
actual stresses applied to the boundary ; whence the conditions to be satisfied 
along a free edge parallel to the axis of y are 

dH x 



N* 



dy 



0, 



and along one parallel to the axis of x, 



N v 



dH x _ 

dx 



0, 



(16) 



(17) 



Since the frequency of the flexural vibrations depends upon the thickness, 
the terms on the right-hand sides of the last two of (10, A) may be omitted in 
the values of N x , N t furnished by these equations, and (16) becomes 

d J d*w 3m — n d?w ) n 

dx \ dx % 2m dy % ) ~~ 
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When the plate is bounded by any plane curve, the two boundary condi- 
tions may be obtained in a similar manner. The work is given on pp. 42 and 43 
of my paper on plane plates (Proc. Lond. Math. Soc, Vol. XXI), and will be 
found to furnish the same results as those given by Lord Rayleigh (Theory of 
Sound, Vol. I, p. 297). 

The preceding analysis gives the complete theory of the flexural vibrations 
of a plane plate. 

JExtensional Vibrations. 

11. We must now consider the extensional vibrations of the plate. The 
value of T x is 

*=£{P+K(% )+**■(£)} «=2 & P +P .(« > ) (20) 

Now P = (m + n) a x + (m — n)(a s + <r 3 ) . 

Substituting the value of a 3 from the first of (5) we obtain 

P=2n{a 1 + E(a 1 + a 2 )\ + EA . (21) 

We therefore see that T x contains the term 2EA Q h, and since by the funda- 
mental hypothesis A is (or may be) proportional to h % , the value of T x cannot 
be obtained by direct integration unless the value of A is known. 

The value of T x can, however, be obtained by an indirect method. I have 
shown in my paper on Plane Plates* that the various constituents of the equa- 
tion which expresses the Principle of Virtual Work can be completely calcu- 
lated as far as the terms proportional to h 3 without knowing the values of the 
quantities A Qt A 1 , A % or the values of o a and a 2 ; and by working out this equa- 
tion by the ordinary methods of the Calculus of Variations the correct values 
of T lf T% and also M can be obtained as far as the term involving h 3 . The 
problem is therefore capable of being completely solved to the order of 
approximation adopted, but for the details of the analysis the reader is referred 
to my paper. In the present article I shall only consider those portions of the 
extensional terms which are proportional to the first power of the thickness. 

Let & = <J 1 + E{p 1 + e t ), % = a z + E(<s 1 +c i ), (22) 

where <r t , a z denote the values of the extensions of a point on the middle surface ; 

*Proc Lond. Math. Soc, Vol. XXI, p. 83. 
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also let a denote the value of a' s at the middle surface. Then from (20), (21) 
and (22) it follows that the value of the term in T x which is proportional to h is 
4nh^, ; and by a similar process we can find the values of T % and if to a first 
approximation. We thus obtain 

T 1 = 4mK&, T 2 = 4nK§, M=2nfa, 

where ^ and |f| are given by (22) and 

du dv dv du /„„.. 

<Tl = ^' G * = -dJ> 6,= ^" , -%- (23) 

Accordingly the first two of (l) become 



(24) 



\dx ^^ dyJ' 
o(v-Y) = 2n^ + i ^). 

The boundary conditions at a free edge are that T x , T 3 and M should vanish 
there. 

The motion which is represented by these equations consists of a displace- 
ment parallel to the middle surface, but there is no displacement perpendicular 
to this surface. It also follows from (23) that the motion involves extension 
of the middle surface, but it does not involve flexion, for the terms upon which 
flexion depends are zero, when w, the normal displacement, is zero. 

From the preceding analysis it follows that the vibrations which- a plane 
plate is capable of executing consist of two independent types, viz. flexural 
vibrations unaccompanied by extension, and extensional ones unaccompanied 
by flexion. 

12. The reader who compares §§7 to 11 with Chapters XIX and XX of Mr. 
Love's treatise cannot, I think, fail to be struck with the simplicity and elegance 
of the above method ; moreover, the complete theory of the extensional vibra- 
tions as far as the terms involving h 3 can also be obtained by employing the 
Calculus of Variations, which Mr. Love's method cannot apparently effect. At 
any rate he has not attempted any such theory in his book. 

Energy of the Plate. 

13. I have shown on page 46 of my paper on Plane Plates that the potential 
energy consists of three distinct sets of terms. The first term is a quadratic 
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function of the extensional strains ct lt a 2 , a multiplied by the thickness. The 
second term is a similar function of the terms on which flexion depends, multi- 
plied by the cube of the thickness. The third term is a quadratic function of 
the extensional strains and their differential coefficients, which is also multiplied 
by the cube of the thickness. To obtain the complete expression for the potential 
energy involves nearly three octavo pages of work which I do not propose to 
reproduce, as the reader who is curious upon the point can consult my paper, 
but the first two terms can be easily obtained by the preceding methods. 

If W be the potential energy of the element 2hdS, we have by the ordinary 
formula 

W = 4 r \{m + n)(a{ + a' t + a' 3 f + n [of + 5J' + of- 4 (<tfrj + oft + *&)] \ dh '- 

By the fundamental hypothesis the terms involving a{, a' % are proportional 
to A 5 and are therefore to be neglected.' The term involving h, upon which the 
purely extensional vibrations depend, can be easily shown to be 

2nh\a\ + al + E(a 1 + cr,) 2 + 1« 2 }, 

where <r 1( o. 2 , a are given by (23), whilst the term in h s which depends upon w is 

Although the term in the potential energy which depends upon bending is 
multiplied by h s , whilst that which depends upon extension is multiplied by h, 
the former term is the most important. One reason of this is, that the gravest 
notes yielded by a plate depend upon flexural vibrations, whilst those depending 
on extension are of high pitch. Also, it is far easier to bend a thin plate than 
to stretch it ; consequently, if any system of forces were applied to the edges of 
the plate, the work done by the forces would be almost entirely concentrated in 
a deformation involving flexion. The extensional part of the deformation would 
be hardly appreciable in the case of a glass or metal plate ; it would only become 
important in the case of highly extensible substances, such as india-rubber. 

Cylindrical Shells.* 

14. The deformation of a naturally curved plate, which will be called a 
shell, is a problem of greater difficulty, because the vibrations are not in general 
separable into two distinct types one of which depends upon the normal dis- 

* See Phil. Trans., 1890, p. 438. 
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placement whilst the other depends upon the tangential displacement. The 
quantities which define the state of strain consist as before of two distinct sets 
of terms, viz. quantities upon which the extension of the middle surface depends, 
and quantities upon which the bending depends. The potential energy, unlike 
that of a plane plate, consists of four distinct sets of terms,* viz. (i) a term pro- 
portional to the thickness multiplied by a quadratic function of the extensions 
of the middle surface ; (ii) a term proportional to the cube of the thickness 
multiplied by a quadratic function of quantities upon which pure flexion 
depends ; (iii) a term proportional to the cube of the thickness multiplied by the 
product of the flexural and extensional terms ; (iv) a term proportional to the 
cube of the thickness multiplied by a quadratic function of quantities upon 
which extension depends. The correct expression for the potential energy of a 
circular cylinder is given by equation (24), page 443, of my paper in the Philo- 
sophical Transactions for 1890, and the values of the various quantities involved 
are given on pp. 440 and 441. 

15. In the figure to §2, let OA be the generator at any point of a right 
circular cylindrical shell, and let OB be the circular section through 0, whilst 
00 is the normal. Let u, v, whe the displacements of along OA, OB and 00 
respectively. Let the unaccented letters denote the values of the different quan- 
tities at 0, and the accented letters at a point J" of which is the projection. 

We shall commence by calculating the values of the four couples. The six 
equations connecting the strains and displacements are 

du' 
dz 



°i = 



. _ 1 /do' , A 



, dw 1 



wi = 



dv> 



dr 



-, dw 



•d 1 dw' 

r ' r d<p ' 
du' 



(1) 



a' 3 = 



+ dr' 
du' dv' 
dty dz 



where r, <£>, z are cylindrical coordinates. 



*Phil. Trans., 1890, p. 443, equation (24). 
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Let A, (i, p denote the values of da 1 /dr, do^/dr, da 3 /dr when r = a. 
The value of Q x is 

G 1 = — f h Q>h'dh' 



(2) 



whilst 



\drJ '. 
Q t =fp (l + -^-) AW 

~1(£ * 



=»4r^) + 4}. 



(s) 



Now 



i? = ( TO + ») ff { + ( m - n )(a{ + <$ 
= (m + n) <y g + (rn — ^(cTi + <T g ) 

+ | (m + n)(^) + (m - n) (A + p) } # + . . . . , 

also iE' = A + ^' + 

where by the fundamental hypothesis A , A x are proportional to h* and h 
respectively ; whence 



cr.= 



A n 



rn -+- n 
/da 3 \ _ A x 
\dr ) rn -\- n 

The values of P and 0/ are 



E{2, + (i). 



(4) 



P = (m + n) o[ + (m- n)(a' t + <r 3 '), 
Q = ( m + n ) <ri + (w - »)(otf + <r{) . 

Expand these values of J", <2' in powers of h', and substitute the values of o s and 
(da 3 /dr) from (4); it will then be found that when the resulting values of P, 
(dP/dr) and (dQ/dr) are substituted in (2) and (3), the terms involving A and A t 
are proportional to some higher power of h than h 3 , and may therefore be 
neglected. We thus obtain 



G 1 = — ihtf\p + E(A + p)\, 

G, = inh s \% + E{% + p) + [a, + E( Gl + <r,)]/a}. 



(5) 
(6) 
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Also by a similar process it can be shown that 

H 1 =-inh*(p + ~f), (7) 

5- 2 = lnh 3 p. (8) 

16. The quantities o lt a % , a s are determined by the equations 

„ _ du 1 /dv , \ _ 1 du , dv ,-,,* 

(Tl -W' <T8== ^(# +M 'J' 6,3= V # + w (10) 



and we must now calculate the values of %, u, p. 
From the first and fifth of (1) we obtain 

/ dOi\ da2 dhn 

^-\dr)~lk~ dz* * 

Now by the fundamental hypothesis a % is proportional to h % , accordingly the term 

da 2 /dz when multiplied by h? is proportional to W and may therefore be 

neglected. We thus obtain 

d?w /.,,\ 

*=— a?" (11) 

Similarly from the second and fourth of (1) we obtain 



'&w . \ e 

Lastly, 



1 /d?w \ E . . /10 , 



/do 3 \ __ 1 / d*u \ / d?v \ 1 ^m 
^ — V ^J — "a V^#j + \drdz) ~~ a* ^ ' 

Substituting the values of (dv/dr), (du/dr) from the fourth and fifth of 
(1), and recollecting that the terms in %, a 2 are to be neglected when multiplied 
by h s , we obtain 

P~ a~ dzchp + a dz ~ a? d$ ' 

Equations (11), (12) and (13) combined with (5), (6), (7) and (8) completely 
determine the couples to the above order of approximation. 

It is very important to notice that the complete values of the couples con- 
tain certain terms depending on the extension of the middle surface, as well as 
terms depending on the bending. 
36 



272 Basset: On the Deformation of Thin Elastic Plates and Shells. 

Inextensibh Deformations . 

17. In a large and important number of problems, the extension of the 
middle surface may be neglected. In this case a 1 , <r 2 , a 3 are zero, and the dis- 
placements are connected together by the three equations which are obtained 
by equating the right-hand sides of (10) to zero ; accordingly, since du/dz = 0, 
we have 

«.= £<•+*>($ + -). <«> 

G *=-iarKw + v 0' (15) 

„ TT 4nh 3 / <ffw dv\ 

H ^~ H * = l^KdWp-lbJ- w 

18. We have now obtained the materials which we require for solving the 
problem of the flexural vibrations of an infinitely long cylindrical shell which is 
bent about a generator. 

The equations of motion of the cylinder are 



1 dT , N 

— ill = vv, 

a dq> a 

1 dN T 
a d<p a 

a d$ 



(17) 



where v = 2ha, a being the density. Also the condition of inextensibility is 

dv , 

dj + w = - ( 18 ) 

Substituting the value of Gc from (14) in the last of (17), we obtain the 
value of N\ and if we substitute this value of iVin the first two of (17) and then 
eliminate T, we shall obtain the equation 

4mnh 2 / d 3 d \ a , d?v , x 

^(ffl + nj^+fj^^-^ ' ( 19 ) 

whence putting 

v = Ae ipt+ "*, 
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we obtain 

V ~ 3pa 4 (»i + w)(s iS + 1) ' 



_ 4mnhV(s z -lf (20) 



which determines the period of the vibrations. 

19. Equation (20) is discussed in Lord Rayleigh's Theory of Sound, Vol. I, 
Chap. X, and also by Prof. H. Lamb (Proc. Lond. Math. Soc, Vol. XIX, p. 365). 
When the cylinder is complete, s is any integer, unity excluded ; but if the cross- 
section of the cylinder consists of a circular arc of length 2aa , s will not be an 
integer. In the latter case the values of s must be determined as follows: 
Along the two generators which form the free edges, T, N and G must vanish ; 
in other words, if the middle point of the arc be taken as the origin of q>, these 
stresses must vanish when $ = ±a. The complete solution of (19) involves 
six arbitrary constants, and the boundary conditions give six equations which 
will enable the six constants to be eliminated ; and the resulting equation com- 
bined with (20) will enable s to be eliminated and an equation to be obtained 
which will determine p in terms of known quantities. 

20. The preceding solution raises a curious point which excited a good deal 
of controversy in England about five years ago. Although the displacement 
parallel to the generators is expressly assumed to be zero, yet in this case there 
is a flexural couple Q % about each circular section whose value by virtue of (15) is 

*=-£*($+•)■ < 21 > 

the extensibility being neglected. 

Let us now suppose that the cylinder is of finite length 2z ; we have in this 
case to satisfy the boundary conditions at the two circular edges, and one of 
these conditions is that 6r 2 = , because the flexural couple about a circular 
section must vanish at a free edge. This requires that 

d?w , 

aW z + W = ' 
whence at the edge 

w = A cos $ + B sin $ 
for all values of $. 

Now it is impossible to satisfy this condition for the most general type of 
non-extensional vibrations, and hence the solution fails ; in other words, when 
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the cylinder is of finite length, purely flexural vibrations unaccompanied by 
extension of the middle surface are impossible. 

This was the controversy raised by Mr. Love in 1888 ; but on examining 
the question it appeared to me that it might be possible for some extension to 
take place which was confined to points in the immediate neighborhood of the 
free edge, and which was inappreciable at points whose distance from the free 
edge was large in comparison with the thickness. But in order to satisfactorily 
deal with the question I found it necessary to obtain the correct values of all the 
quantities involved as far as terms proportional to the cube of the thickness, on 
the supposition that the middle surface undenoent extension. This required a 
long and difficult investigation by means of the Calculus of Variations, which is 
given at length in my paper (Phil. Trans., 1890, p. 433) ; but I was finally 
enabled to establish the correctness of my supposition by solving the following 
problem. 

21. Let a very long, heavy cylindrical shell, whose cross-section is a semi- 
circle, be suspended by vertical bands attached to its straight edges, so that its 
axis is horizontal. It is required to investigate the state of strain produced by 
its own weight. 

I found that if R denote the change of curvature along a circular section, so 
that 



„ _ 1 /d 2 w \ 



w being the displacement along the radius, and <r 9 is the extension of the middle 
surface along the circular section, so that 



1 /^ ■ \ 



.d(p 
then 

B E 3a (^ it — 4* sm 4> — cos 4*) 



et % a ~h % \\n — cos<?> + §E($7t — $ sin<£> — cos<£>) [ ' 

where $ is measured from the lowest point. 

Since the numerator of this fraction is an even function. of <£>, it does not 
change sign with q> ; also the numerator is always positive between the limits 
^71 and — in, and its maximum value occurs when $ = 0, and is equal to 



= 


E 

a 


+ 


3a 


B 


— 


E 
a 


• 
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\it — 1 , and its minimum value occurs when $ = \n and is equal to zero. We 
therefore see that when ty =■ , 



and when.^) = \it, 



Since the thickness of the shell is supposed to be small compared with its 
radius, it follows that the change of curvature is large compared with the exten- 
sion of the middle surface except when a(j7t — $) is comparable with h, i. e. in 
the neighborhood of the straight edges of the shell ; and therefore at all points of 
the shell whose distances from the edges are large in comparison with the thick- 
ness, the terms depending on the product of the change of curvature and the 
cube of the thickness, i. e. the terms upon which the bending depends, are of the 
same order as those depending upon the product of the middle surface and the 
thickness ; but at points whose distances from the edge are comparable with the 
thickness of the shell, the extension of the middle surface is of the same order as 
the change of curvature, and therefore the terms depending upon the product of the 
change of curvature and the cube of the thickness are small in comparison with 
those depending upon the product of the extension and the thickness. 

22. Now we have already pointed out in §14 that the potential energy con- 
sists of four distinct terms, the first, third and fourth of which involve extension, 
whilst the second depends upon pure flexion unaccompanied by extension ; and 
the preceding result shows that the second term is by far the most important. 

Under these circumstances we are justified in concluding that the solution 
given in §18 is substantially correct, and that any extension which takes place 
is practically confined to points in the neighborhood of a free edge, and is practi- 
cally negligible. 

23. By means of the fundamental hypothesis that the three stresses B', S', T 
may be treated as zero, the value of the potential energy when the middle 
surface undergoes no extension can easily be shown to be 

W=inh 3 {2?+{i*+E(X + (if + \tf }. 
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Substituting the values of /z, p from (12) and (13) and taking into account 
the conditions of inextensibility — see equation (10) — which require that % = 0, 
this may be expressed in the form 

W=^\ - f^ + wY+f— -— Yl. (22) 

W 3a? \ {m + n) a 2 \dty z ^ w )^ \dzdty dz ) ] ' K ' 

From the second of (10) this may be expressed in terms of v if desired. 
The conditions of inextensibility (10) require that 

u=Ua,v = -z^ + V,w = z^f-^-, (23) 

and since ZJand V must necessarily be periodic functions of ty, provided the 
cylinder is complete, we may express them in series of the form 

where A and B are complex functions of the time. The potential energy of a 
cylindrical shell of length 2z may be obtained by substituting the above value of 
u, v, w, U, V in (22) and integrating throughout its length. Writing out the 
full values of u, v, w in real quantities, we get, as shown by Lord Rayleigh,* 

u ■■=■ — s^al, (B s sin sty + B', cos sty) , 

v — 2-j (A s a + B s z) cos sty — (A' s a + B'^z) sin sty), 

w = 2 ■{ s (A,a + B g z) sin sty -\- s (A' s a + B' s z) cos sty \ . 

The value of the potential energy can now be obtained in terms of the A's 
and B's, whilst the kinetic energy is obtained from the equation 



T— tyffW + & + w a ) adtydz. 



Having evaluated T and W, the motion can be obtained by Lagrange's 
equations, and it will be found that the vibrations consist of two distinct types. 
The first type depends upon the J.'s, and are those which we have already inves- 
tigated ; they involve pure flexion, and are such as would be produced if every 
circular section were slightly deformed into an equal and similar curve. The 
second type involve torsion as well as flexion, and each circular section is twisted 
into a sinuous curve about its mean position. 

*Proc. Boy. Soc, Vol. XLV, p. 105. 
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The periods of the latter class of vibrations are given by the equation 

„_ 4mwftV(g 2 — If l + 3o?(w + n)/ffM , P 
P ~ 3pa 4 (m + n)(«» + 1) ' 1 + 3a 2 /s 2 (s* + 1) P ' 

where s = 2 , 3 , . . . . 

For the further discussion of these vibrations we must refer to Lord Ray- 
leigh's paper ; we notice that if the cylinder be at all long in proportion to its 
diameter, the difference between the periods of the two types of vibrations 
becomes small. 

24. Before leaving this branch of the subject, it may be well to point out 
that Kirchhoff's four boundary conditions apply to the case of a cylindrical shell 
as well as to a plane plate.* The principles upon which the theory rests are 
identical in both cases, although the mathematical details are different. 

Along a free edge which is a circular section, the conditions are 

21 = 0, 6? 2 =0, 

M^HJa, * = 4-tf, 
whilst along a generator the conditions are 

T t —0, M 1 = 0, 

ffl = 0, * 1 + <f =0. 
The proof depends upon the fact that the integral 

when taken round a rectangle bounded by four lines of curvature, vanishes by 
Stokes' theorem. 

Extmsional Vibrations. 

25. The potential energy contains a term proportional to the thickness 
multiplied by a quadratic function of the extensions. The value of this term is 

W= 2nh\af + <x» + j0 ( tfl + <r 2 ) 2 + i*5f , 

which can easily be established by the methods of the present paper. 

*See p. 453 of my paper, Phil. Trans., 1880. 
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The values of T lt T % , M x and M s , so far as they depend on h alone and not 

0n **' are T x = 4nh& , M x - 2nha 3 , 

T % = 4nh^, M 2 = 2nha 3> 
where ^ = a 1 + E{a x +a t ), % = o i + E(c x + a z ) . 

The equations of motion are 

ow = — 2n%/a + pZ, 
whilst the boundary conditions at a free edge are that T lt T 2 , M lt M % should 
vanish there. 

The solution of these equations when no bodily forces act has been dis- 
cussed by Mr. Love,* and also by Lord Rayleigh,| to whose papers the reader 
is referred. The vibrations are of high pitch, and are of mathematical rather 
than physical interest. 

26. Before passing on, it may be as well to take this opportunity of reply- 
ing to a criticism which Mr. Love has made in the notes to pp. 238 and 262 
of Volume II of his Treatise on Elasticity. He observes that the value of T % 
given by the second of my equations (44) (Phil. Trans., 1890, p. 450) does not 
furnish the same value of T % as the equations given at the bottom of p. 456, and 
not being able to understand this, he imagines that he has detected some error 
in my work. The explanation is simple. In the investigation by which equa- 
tions (44) were obtained, the middle surface is expressly supposed to undergo 
extension, whereas in the problem discussed on p. 456 no extension is supposed 
to take place. The two cases require totally different treatment, and are there- 
fore not parallel ; consequently Mr. Love's so-called test is nugatory. If there 
were any error in the value of T % , the equations obtained by substituting the 
values of the stresses given by (43) and (44) in the first three of (11) would not 
reproduce equations (45), (46) and (47), as they actually do. 

The point where Mr. Love has gone astray is the following : 
Let .Fbe any function of u, v and their differential coefficients with respect 
to x and y. The condition that the integral 

CfoFdxdy 

* Phil. Trans., 1888, pp. 538, etc. t Proc. Roy. Soc, Vol. XLV, p. 443. 
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should vanish is obtained by integrating it by parts until it is reduced to the 
form 

/ / (M8u + NSv) dxdy + a line integral. 

Now if Su, $v are independent quantities, the surface integral gives rise to the 

equations 

M=0, N=0, 

which must hold good at every point of the surface, and the line integral gives 
rise to similar equations which must hold good at the boundary. If, however, 
Su, Sv are not independent but are connected together by some relation, as is 
the case when the surface undergoes no extension, the equations M = 0, N= 0, 
as well as those furnished by the line integral, are no longer true. When the 
relation is of the form 

v =/(«, du/dx, ....)> 

so that v is given explicitly in terms of u, we can substitute for hv in terms 

of Su and integrate the term / / Nhvdxdy until it is reduced to the form 

N'hvdxdy -f- a line integral, 



ff> 



and since 8u is independent, the surface integral furnishes the single equation 

M+N' = 0, 

and the equations furnished by the line integral terms will also be different. 

Instead of adopting the above-mentioned procedure, we may, if we please, 

introduce an indeterminate multiplier, in which case we shall obtain the two 

equations 

M + p = 0, N+q=0, 

where if and N have their former meanings and p and q are additional terms 
which depend upon the indeterminate multiplier. The elimination of this 
quantity must necessarily lead to an equation equivalent to M + N 1 = 0. 

The failure to recognize the distinction between the two cases constitutes 
the error of Mr. Love's criticism. 

Spherical Shells.* 

27. The theory of spherical shells can be worked out in exactly the same 
manner. In the figure to §2, let OA be a meridian and OB a parallel of lati- 

*Phil. Trans. 1890, p. 463. 
37 
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tude. The six strains are given by the equations 

, 1 / du' 



f^W + W 0' 



— ( - — 2 ~j^r + w' cot + w' ) , 
r \sm d<p ' J 



oi 



5{: 



dud 
dr 

1 dvd dv? 
r sin 6 d<p dr 



r 



dvl 



u 



1 dw' 



a *~ dr~ r + r d$ ' 
_ __ 1 { del 
"s-TKdd 



Vcot0-j- 



du' 



\ 



sin 6 d$ /' 
By direct integration we find 






-*•*{(§)+-!}• 



= *, + <# fa + <t 8 ), 
= ^ + #(*, + <«) . 



(1) 



(2) 

(3) 
(4) 



Hi — — s % 

Putting as before 

& = o 1 + E(<y 1 + o< i ), 
%=*, + E(% +[i) , 
equations (2), (3) and (4) become 

<* = -♦»*■(£ + A), G 2 =$nh>(%+^-),) 

where %, fi, p denote the values of ddijdr, do % jdr and da s /dr when r — a. 

Recollecting that B', o{ and 51 are to be treated as zero by virtue of the fun- 
damental hypothesis, we readily obtain to the required order of approximation, 
1 / d?w \ E 



(5) 
(6) 



(7) 



1/1 d?w , . a dvo . \ E . 

f = ~^K^Qd^ + COte W jrW )—a-^ + a " ) ' 
2 / , g dw d?w \ 
a 2 sin V <fy ~* dfcty/ ' 



3,= 



jp = 



(8) 
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The proof of the above formulae may be left to the reader, and he will find 
no difficulty in establishing them by means of the methods explained in the earlier 
part of this paper. 

28. We shall now calculate the potential energy in the particular case in 
which the middle surface undergoes no extension. 
The expression for Wis 

W = \ffj\m + n) A' + n{ a" - 4 (a'A + ofa + aft) }] (1 + h'/dfdh'dS, (9) 

where A' is the dilatation at the point a + h 1 , 0, ^> ; a' is written for a' 3 , and 
o'i, of are omitted because they would give rise to terms proportional to h*. 

Since there is no extension, A — A /(m-\-ri), and since A is proportional 
to A a , the only term in A' which it is necessary to retain is h! (dA/dr) ; we accord- 
ingly obtain 

i (m + n)J_A>° (1 + h'/afdh' = ~^ (X + tf ; 

also since o^ = a 2 = , we have 

a[ = h'%, a^h'fi, 
so that Znfolol (l + h'/af dh< == \ fttfty . 

In the same manner 

2n fjp{ + a'*) ai (1 + h'/af dh = —% nh 3 (X + (if 
afld i nfl" (1 + h\ of dh! = \ nh 3 p*. 

Hence if TFnow denote the potential energy per unit of area, 

W=inh s {tf + (i* + E(h + ny+lp>\. (10) 

The values of X, {i, p are given by (8), but the extensional terms must be 
put equal to zero. 

The complete expression for the potential energy as far as A 3 , when the 
middle surface is supposed to undergo extension, is given by equation (16), 
p. 467 of my paper in the Phil. Trans., 1890. It will be observed that this 
expression (as in the case of a cylindrical shell) consists of four distinct sets of 
terms : (i) a term proportional to h multiplied by a quadratic function of the 
extensions, whose value is 

2nh\<r\ + a* + E(o l + o,Y + %(>*}; (11) 
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(ii) a term proportional to h 3 , which is multiplied by a quadratic function of 
quantities upon which pure bending depends ; (iii) a term proportional of h 3 
multiplied by the products of extensional and flexural terms ; (iv) a term pro- 
portional to h s multiplied by a quadratic function of extensional terms. It will 
not, however, be necessary to consider the last two terms. 

Flexural Vibrations. 

29. Let p x , p 8 be the radii of curvature of the deformed surface along and 
perpendicular to a meridian ; then 



1 /dfw , \ 



JL J_ _ JL ( dS>w 

Pi ~~ a ~ ~~ 

1 1 1/1 d % w , .A^tV \ 



(12) 



(13) 



p 2 a 

These formulae are probably well known in America, but if a proof is 
desired one will be found on p. 477 of my paper. "We notice from (8) that when 
extension does not take place, X and ft are proportional to the changes of curva- 
ture along and perpendicular to a meridian. 

The flexural vibrations of a spherical bowl have been investigated by Lord 
Eayleigh.* The conditions of inextensibility are given by the first, second and 
sixth of (1) and are 

du , 

-w + w==0 > 

dv , „ , ■ a n 

-rr- + U COS + W Sin = , 

du a , ■ a dv 

^-<»cos0+sm0 w = O.J 

Let % = log tan \ 0, then since 

. n d d 

we find upon eliminating w that 

d / u \ d / v \ 

d$ Vsrnlj + di V sin 0/ ' 

^Vsin0y ^Vsin0y"~ u ' 

consequently u cosec and v cosec both satisfy an equation of the form 

dTV ,dTV_ 
# 2 + d z 2 -°- 



*Proc. Lond. Math. Soc, Vol. XIII, p. 4. 
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In the case of a spherical bowl, all the quantities are necessarily periodic 
functions of <£>, and we may therefore suppose that $ enters in the form of the 
factor e'** where s is an integer ; whence 

dtf ' 

giving F=^tan 8 i0 + £cot s £0. 

If is measured from the pole, we must have B = , otherwise V would 
become infinite when 0=0; we may therefore take 

u= — 24 s e' s *sin0tanH0, \ 

v= t2.4,e' 8 *sin0tan s £0, I (14) 

w = 2 A s s 1 '* (s + cos 0) tan* \ , ) 

where 8=2,3,4,...., and A, is a complex function of the time. From these 
equations it can be shown by (12) that 

JL__L = _ s ±£ = 4£Mi! = _(J___L') > (15) 

also by the last of (8) , 

_ n ^ A,(s s -s)s>°*ta,n>$6 , . 

P=~^ a-^B * (16) 

The value of A, may be supposed to be real without loss of generality ; 
accordingly, rejecting the imaginary part, we obtain from (10), (12), (15) and (16), 

W = ®?J- 3 /"cosec 3 $X Ay (s—lf tan 8 * \ dd , 
oa •'o 

where a is the angle of the bowl. 

Lord Kayleigh has applied the above expression for W to determine the 
notes produced by a hemispherical bell, but he has not attempted to take into 
account the conditions which must be satisfied at a free edge. Now one of these 
conditions is that the flexural couple Q % should vanish there, where 

G i =inh 3 \^ + E(^ + l i)\ 

V pi a / 

From (15) we see that O z cannot vanish at a free edge, and consequently 
a spherical bell whose edge is free cannot vibrate in this manner if the middle 
surface is supposed to remain absolutely inextensible during the motion ; but 
since the complete value of G% contains terms depending on the extensibility, it 
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would be possible to satisfy the boundary conditions if extension were supposed 
to take place. From the corresponding results in the case of a cylindrical shell, 
we should expect that some local extension must necessarily take place in the 
immediate neighborhood of a free edge ; but at points whose distance from this 
edge is at all large compared with the thickness the extension becomes negligible, 
so that the solution given by Lord Rayleigh is for all practical purposes a suffi- 
ciently accurate one. 

Extensional Vibrations. 

30. We have already pointed out that the potential energy contains a term 
proportional to the thickness which depends upon the extension alone. This 
term leads to the following equations, viz. 

?i=4wAJt, T i = inh'§, M 1 = M z =2n7ia, (17) 

and the equations of motion are 






(18) 



Equations (17), which are approximate, may be deduced by direct integra- 
tion, whilst (18; may be obtained by writing down the equations giving the 
motion of translation of an element. The terms involving the normal shearing 
stresses N lt N 2 must be omitted because these quantities are proportional to the 
cube of the thickness. Having done this, we substitute the values of T lt T % 
and M from (17). 

The integration of these equations has been discussed by Mr. Love. 

Collapse of a Boiler Flue. 

31. A problem of great practical importance in engineering is that of the 
collapse of a boiler flue. The pressure of the smoke and the hot air from the 
furnace which passes through the flue is approximately equal to the atmospheric 
pressure, whilst the pressure of the steam within the boiler is considerably 
greater. Under these circumstances the flue shows a tendency to collapse when 
the difference between the pressures is sufficiently great. 
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The complete mathematical theory has not yet been made out, and the diffi- 
culty arises from the fact, that since the flue is subjected to an external pressure 
7ii, where % is the pressure of the steam in the boiler, and an internal pressure 
7t 2 , where tl% is the pressure of the hot air from the furnace, it is not permis- 
sible to suppose that B' may be treated as zero throughout the substance of the 
flue ; for as we pass from the outer to the inner surface B' must vary from 
— 7ti to — n z . The value of B' is 

B' = A + AJt! + i A 2 h'*+ . 

where h' is a small quantity, but the A's are not ; we must therefore retain 
them in the expression for the flexural couple, and consequently the value of this 
quantity cannot be found unless the A's are known. 

In a paper published in the Phil. Mag. for Sept. 1892 I fully discussed the 
difficulties of constructing a satisfactory theory. It seems allowable to assume 
that S and T may be treated as zero, in which case the value of the flexural 
couple G about a generating line can be shown to be 

(x -~ f£lAl/l ~5a(m + n)\m + n~ a ~ a dp J' 

the extension being neglected. 

As we do not know anything about the A's, I assumed provisionally that 

n W / rnn \/d*w \ 

where a is some unknown function of the pressures n 1 and 7%, which vanishes 
when these quantities are zero. 

The equations determining the small oscillations are 

dT N ol 

ds ■ p ' 

dG , AT n 

where p is the radius of curvature of the deformed middle surface and 0" the 
density. The condition of inextensibility is 

£ + .-0, (2) 
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also 



o - a ~ a* W / 



p a a* \d<p* 
and ds = adty . The equations of motion may now be written 

~^+N=2haav, 

dN m . , s \ 7 , , X /^ 2 w 






Substituting the value of G from (1), eliminating T and JV and taking 
account of (2), we finally obtain 

{la- 1 (Z> 2 + 1) + Th — ih]^ + 1) D% + 2ha<r (Z? 2 — 1) y = 0, (3) 

'=*?(=?= + •)• " = <"*• 

To solve (3) assume that 

where s = 2, 3, 4, then 

\Ia~ l (s z — 1) - (Ttj - i%)\(s 2 — 1)«» = 2£a<r (s 8 + l)p>, 

and therefore j? will be real provided 

7ti — Tt% < /(s 8 — l)/a. 
The least value of the right-hand side occurs when s = 2, in which case 

7ti — ?t a < 3//a 

^ 8h 3 / mn \ , ,. 

<-r( ■ +a). (4) 

If it were allowable to neglect a, the condition would become 

s Smnh 3 /.s 

a 3 (m + n) 

Finite Bending of Thin Shells. 

32. The problem of the finite bending of a thin shell presents difficulties 
which cannot yet be said to have been completely surmounted. When the 
deformation is small, the displacements of a point on the middle surface are 
small quantities whose squares and higher powers may be neglected in the cal- 
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dilation of the strains ; but when the deformation is finite the displacements 
may be large quantities. Since any deformation involves a change in the values 
of certain geometrical quantities, such as the curvature and torsion of certain 
lines drawn on the middle surface, the most appropriate course to pursue is to 
endeavor to express the stresses in terms of such geometrical quantities. 

There is one class of problems which can often be solved without much 
difficulty which arise when a plane surface is bent without extension into a 
developable surface, or when a developable surface is bent into a plane or into 
some other developable surface such that the lines of curvature on the old surface 
are lines of curvature on the deformed one. This method can generally be 
applied when a plane plate is bent into a conical or cylindrical surface.; but it 
could not be applied when a right circular cone is bent into a cone whose lines 
of curvature are not identical with those of the original one. 

The success of this method, in cases where it can be applied, depends upon 
the fact that the flexural couples Gr lt Gr % can be expressed in terms of the 
changes of c rvature, and also that in the special cases alluded to a sufficient 
number of the ten stresses are zero to enable the remainder to be determined by 
means of the general equations of equilibrium. 

33. Let OA, OB be two lines of curvature on the middle surface of the 
undeformed shell; lt 2 the centres of principal curvature; let oa, oh be the 
curves in which the planes OA O x , OB0 2 meet any layer of the shell. Let 






p 1; p 2 be the principal radii of curvature at ; let 0o = ?/; also let accented 
letters denote the strained positions of the various points. 

If P denote the normal traction along oa, and B the normal traction 
along Oo, P = (m + n)a 1 + (m — n)(a % + a,) 

= 2n(<7 1 + E{g 1 + a,) \+EB. (1) 

38 
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Now _ o'a' — oa 

a±— 

oa 
and oa . , V o'a' yf . , . 

Since we neglect the extension of the middle surface, O'A' — OA, whence 

=<i-i) + ^{^-^ + *4 (2) 

Similarly, 

and the value of G 2 is 

G t = f h Pridri. (4) 

Now, according to the fundamental hypothesis laid down in §5, it follows 
that, provided there is no external pressure, B must be a quadratic function of h 
and r\, and consequently the retention of B will lead on integration to terms in 
Q % involving higher powers of h than h 3 ; accordingly we may neglect B, and if 
we substitute the values of <r 1( a z from (2) and (3) in (1), and the resulting 
value of P in (4) and integrate, we shall obtain 

Similarly we can show that 

*=-*•*{*(■£-£) +< i+ *>0s--i)}- (e) 

We notice that the above values of Gr x and 6? 2 agree with those which we 
have already obtained for plane plates and cylindrical and spherical shells when 
the bending is small. 

When a plane plate is bent into a developable surface p x = p 3 = co ; also one 
of the quantities p{ or o 2 (say p0 is infinite; whence (5) and (6) become 

G 2 = Uh s (l + E)/p[,\ , 

Since the extension of the middle surface is neglected, equations (5) and (6) 
would not apply to the case of a plane plate which is deformed into a surface 
such as a portion of a sphere. 
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34. As an example of the preceding method, we shall consider the case 
of a plane plate of thickness 2h which is bounded by two radii OA, CD and 
two arcs OB, AD of concentric circles; and we shall inquire whether it is 
possible to bend this plate into a portion of a right circular cone in which OA , 
BD are generators and OB, AD are circular sections. 

"We shall assume for trial that the bending may be effected by means of 
tensions, normal shearing stresses and flexural couples applied to the edges ; so 
that the tangential shearing stresses and the torsional couples are zero. 

A 




whence 



Let a be the semi-vertical angle of the cone, r the distance of any point on 
OADB from 0. Then 

Pi=p2=°°; 9%—^, pi = r tan a; 
G % = inMEr- 1 cot a, 
<?! = — f nh 3 (1 + E) r- 1 cot a. 

The equations of equilibrium are 

dT, 
d<p 

d f%.T ^ • . dN, m 
^- (N 2 r) sin.a + -^ — T % cos a = , 

dG x 



+ N-i cos a = , 



Let 
then 
whence 

and therefore 



£(G a r)-N,r+G 1 = 0. 

inh 3 (l + E)cota = k, 

G^-k/r; 

N t = — k/r», N 1 = 0, 

T 2 = k tan a/r 2 , 



(8) 



T,= 



A k tan 



a 



r* 



where A is the constant of integration. 
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From these equations we see that all the stresses are determinate except 1\. 
If CO = a, OA — b, and T a , % denote the tensions along OB, AD, we have 

A h tan a 



T = 



a a? ' 



T A k tan a 

6_ T~ F~ ' 

from which it appears that either T a or T 6 may if we please be made zero, pro- 
vided the other is properly determined. 

The above results would also apply to a belt of a complete cone bounded 
by two circular sections. 

Fledborough Hall, Holyport, Berks, England. 



[Corrections bt A. Chessin.] 

Page 165. 

Instead of: for every whole negative number a e there should be another 
whole negative number a e ,, and that at the same time A e x a <' + A e ,x ae = 

read : 
a e , a ei , a e being the negative numbers, the expression 

A 6i x~ a H + A e ^r a ^ + + A 6 x~\ 

vanishes identically. 

Page 182. 

Instead of: for every such a e there should be another whole negative num- 
ber a 6 >, and that at the same time 

A e x- a * A e ,x- a °' _ n 



(—aa. + b)l^(— aa 6 , + b)\ 

read : 
o^e,! a e a > • • • . , « e being the negative numbers, the expression 



"p 



A e x- a °> A^x- a '> Ae^xS 

~t~ / _„ L l\ I T " " T 



(— aa ei + b)\ ^ (-aa 6l + b) ! "•" * * ' ' T (- oo. p + 6)! 
vanishes identically. 



